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Abstract 

Hydrodynamic  features  of  a  falling  cylinder  into  the  water  column  are  investigated  experimentally.  The  experiment  consisted  of 
dropping  three  cylinders  of  various  lengths  into  a  pool  where  the  trajectories  were  filmed  from  two  angles.  The  controlled  para¬ 
meters  are  cylinder’s  physical  parameters  (length  to  diameter  ratio,  center  of  mass  location),  and  initial  drop  conditions  (initial  velo¬ 
city,  and  drop  angle).  Six  trajectory  patterns  (straight,  spiral,  flip,  flat,  seesaw,  combination)  are  detected  during  the  experiment.  The 
center  of  mass  position  has  the  largest  influence  on  the  trajectory  of  cylinders.  The  observed  motion  of  cylinders  is  well  simulated 
using  a  numerical  model  based  on  the  momentum  and  moment  of  momentum  balances  using  triple  coordinate  transform. 
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1.  Introduction 

Study  on  the  movement  of  a  rigid  body  in  fluid  has 
wide  scientific  significance  and  technical  application.  It 
involves  nonlinear  dynamics,  flight  theory,  body-fluid 
interaction,  and  instability  theory  (e.g.,  [1]).  The  techni¬ 
cal  application  of  fluid  mechanics  of  a  rigid  body  in  fluid 
includes  aeronautics  and  navigation.  Recently,  the  scien¬ 
tific  problem  about  rigid  body  movement  in  the  water 
column  drew  attention  to  the  naval  research.  This  is 
due  to  the  threat  of  mines  in  the  naval  operations.  With¬ 
in  the  past  15  years  three  US  ships,  the  USS  Samuel  B. 
Roberts  (FFG-58),  Tripoli  (LPH-10)  and  Princeton 
(CG-59)  have  fallen  victim  to  mines.  Total  ship  damage 
was  $125  million  while  the  mines  cost  approximately  $30 
thousand  [2].  Mines  have  evolved  over  the  years  from 
the  dumb  “horned”  contact  mines  that  damaged 
the  Tripoli  and  Roberts  to  ones  that  are  relatively 
sophisticated — non-magnetic  materials,  irregular  shapes, 
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anechoic  coatings,  multiple  sensors  and  ship  count 
routines.  Despite  their  increased  sophistication,  mines 
remain  inexpensive  and  are  relatively  easy  to  manufac¬ 
ture,  keep  and  place.  Water  mines  are  characterized  by 
three  factors:  position  in  water  (bottom,  moored,  rising, 
and  floating),  method  of  delivery  (aircraft,  surface,  sub¬ 
surface)  and  method  of  actuation  (acoustic  and/or  mag¬ 
netic  influence,  pressure,  contact,  controlled). 

Prediction  of  a  falling  rigid  body  in  the  water  column 
is  a  key  component  in  determining  the  impact  speed  and 
direction  of  mine  on  the  sediment  and  in  turn  in  deter¬ 
mining  its  burial  depth  and  orientation.  In  this  study, 
a  cylinder  drop  experiment  was  conducted  to  investigate 
dynamical  characteristics  of  the  falling  cylinder  through 
the  water  column. 


2.  Triple  coordinate  systems 

Consider  an  axially  symmetric  cylinder  with  the  cen¬ 
ter  of  mass  (COM)  X  and  center  of  volume  (COV)  B  on 
the  main  axis  (Fig.  1).  Let  (L,d,x)  represent  the  cylin¬ 
der’s  length,  diameter,  and  the  distance  between  the 
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Nomenclature 

(On,  Cd 2)  drag  coefficients  along  and  across  the 

Re 

Reynolds  number 

cylinder 

y 

translation  velocity  (ms-1) 

(Ctdi,  Ctd2)  translational  drag  coefficients  (kgs-1) 

vr 

water-to-cylinder  velocity  (ms-1) 

C\  lift  coefficient 

Vi 

component  of  Vr  along  the  cylinder  (ms-1) 

Ct  1  translational  lift  coefficient  (kgs-1) 

V2 

component  of  Vr  perpendicular  to  the  cylin¬ 

d  cylinder  diameter  (m) 

der  (ms-1) 

(/ib/2,/3)  added-mass  ratios  for  drag  and  lift  forces 

Vw 

water  velocity  (ms-1) 

/r  added-mass  ratio  for  moment  of  drag  and  lift 

y(m) 

initial  speed  of  dropping  cylinder  (ms-1) 

forces 

V* 

nondimensional  translation  velocity 

(/rd2/rd3)  rotational  drag  force  (N) 

(«o,- 

. . ,  a4)  correlation  coefficients  for  predicting  bot¬ 

Fb  buoyancy  force  (N) 

tom  impact  velocity  and  orientation 

Fd  drag  force  (N) 

[j&o(0,...,j84(0]  correlation  coefficients  for  predict- 

(Fdl,  Fd2,  Fd3)  drag  force  in  the  F-coordinate  (N) 

ing  translation  velocity  and  orientation  of 

Fi  lift  force  (N) 

the  falling  cylinders 

(Fn,  Fi2,  Fi3)  lift  force  in  the  F-coordinate  (N) 

<5 

aspect  ratio  of  the  cylinder 

(}e,  j e,  kg)  unit  vectors  in  the  F-coordinate 

A 

distance  between  adjustable  copper  cylindri¬ 

(ip,  ]F,  kp)  unit  vectors  in  the  F-coordinate 

cal  part  (m) 

(iM,  \M,  kM )  unit  vectors  in  the  M-coordinate 

V 

molecular  viscosity  of  the  water  (m2s-1) 

(Ji,  /2,  J3 )  moments  of  gyration  (kgm2) 

n 

volume  of  the  cylinder  (m3) 

(j\\  4°,  4^)  moments  of  gyration  for  cylindrical 

p 

density  of  the  cylinder  (kgm-3) 

part-z  (kgm2) 

Pw 

density  of  the  water  (kgm-3) 

L  length  of  the  cylinder  (m) 

X 

distance  between  COM  and  COV  (m) 

(, l\ ,  /2,  h)  lengths  of  the  cylindrical  parts  (m) 

Ml 

i/t 2,  1/^3)  angles  determining  the  cylinders’ 

(rn\,..  -  ,m6)  masses  of  cylindrical  parts  (kg) 

orientation 

Mb  torque  due  to  the  buoyancy  force  (kgm2  s-2) 

CO 

angular  velocity  (s-1) 

Mb  torque  due  to  the  hydrodynamic  force 

(a>u 

co2,  co3)  angular  velocity  components  in  the  M- 

(kgm2  s~2) 

coordinate  (s-1) 

(Mdl,  Md2,  Md3)  torques  due  to  the  drag  force  in  the 

K> 

co 2,  (0 3)  angular  velocity  components  in  the  F- 
coordinate  (s-1) 

M-coordinate  (kgm2s-2) 

(Ri,  F2,  F3)  radii  of  cylindrical  parts  (m) 

1 

1  ^ 

: 

S'  X 

1 

-  L/2  - 


-  \J2.  - 


Fig.  1.  M- coordinate  with  the  COM  as  the  origin  X  and  (im,  jm)  as  the 
two  axes.  Here,  %  is  the  distance  between  the  COV(F)  and  COM  ( L,d) 
are  the  cylinder’s  length  and  diameter. 


two  points  (X,B).  The  positive  ^-values  refer  to  nose- 
down  case,  i.e.,  the  point  X  is  lower  than  the  point  B. 
Three  coordinate  systems  are  used  to  model  the  falling 
cylinder  through  the  air,  water,  and  sediment  phases: 
earth-fixed  coordinate  (F-coordinate),  main-axis  follow¬ 
ing  coordinate  (M-coordinate),  and  force  following 
coordinate  (F-coordinate)  systems.  All  the  systems  are 
three-dimensional,  orthogonal,  and  right-handed.  The 
origin  of  M-  and  F-coordinates  is  located  at  X  [3]. 


Fig.  2.  Three  coordinate  systems. 

The  F-coordinate  is  represented  by  FE  (O,  i,  j,  k)  with 
the  origin  ‘0\  and  three  axes:  x-,  y-axis  (horizontal)  with 
the  unit  vectors  (i,j)  and  z-axis  (vertical)  with  the  unit 
vector  k  (upward  positive)  (see  Fig.  2).  The  position  of 
the  cylinder  is  represented  by  the  COM  position 
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X  =  x\  +  y\  +  zk.  ( 1 ) 

The  translation  velocity  is  given  by 

^T  =  V,V=  (u,v,w).  (2) 

The  orientation  of  the  cylinder’s  main-axis  (pointing 
downward)  is  given  by  iM.  The  angle  between  \M  and 
k  is  denoted  by  \j/2  +  n/2.  Projection  of  the  vector  \M 
onto  the  (x,y)  plane  creates  angle  (ij/f)  between  the  pro¬ 
jection  and  the  x-axis  (Fig.  2).  The  M- coordinate  is  rep¬ 
resented  by  Fm  (X,iM,jM,kM)  with  the  origin  ‘X’,  unit 
vectors  and  coordinates  (xM,yM,zM).  The 

unit  vectors  of  the  M-coordinate  system  are  given  by 
(Fig.  2) 

jj»=kx  iM,  kM  =  iM  x  jM  (3) 

The  M-coordinate  system  is  solely  determined  by  the 
orientation  of  the  cylinder’s  main-axis  \M. 

The  F-coordinate  is  represented  by  Ff  (X,ip,  jp^kp) 
with  the  origin  X,  unit  vectors  (iF, jp-,kp),  and  coordi¬ 
nates  (xF,yF,zF).  Let  Vw  be  the  fluid  velocity.  The 
water-to-cylinder  velocity  is  represented  by 

vr  =  Vw  -  V,  (4) 

which  can  be  decomposed  into  two  parts, 

Vr  =  Vi+V2,  Vi  =  (Vr  •  ip)ip, 

V2  =  Vr-(Vr-ip)iF,  (5) 

where  Vi  is  the  component  parallel  to  the  cylinder’s 
main-axis  (i.e.,  along  iM),  and  V2  is  the  component  per¬ 
pendicular  to  the  cylinder’s  main-axial  direction.  The 
unit  vectors  for  the  F-coordinate  are  defined  by  (column 
vectors) 

if  =  *m,  jF  =  V2/|V2|,  kf  =  ifXjF.  (6) 

In  the  F-coordinate,  the  hydrodynamic  forces  (drag,  lift) 
and  torques  are  easily  computed  [4-6]. 


3.  Dynamics 

3.1.  Momentum  balance 


The  translation  velocity  of  the  cylinder  (V)  is  gov¬ 
erned  by  the  momentum  equation  in  the  F-coordinate 
system  [1] 


u 

'o' 

V 

=  - 

0 

w 

_g_ 

Fb+Fh 

Pn 


(7) 


where  g  is  the  gravitational  acceleration;  77  is  the  cylin¬ 
der  volume;  p  is  the  rigid  body  density;  p77  =  m ,  is  the 


cylinder  mass;  Fh  is  the  hydrodynamic  force  (i.e.,  surface 
force  including  drag,  lift,  impact  forces);  Fb  =  - pw77 ,  is 
the  buoyancy  force;  and  pw  is  the  water  density.  The 
drag  and  lift  forces  are  calculated  using  the  drag  and  lift 
laws  with  the  given  water-to-cylinder  velocity  (Vr).  In 
the  F-coordinate,  Vr  is  decomposed  into  along-cylinder 
(Vi)  and  across-cylinder  (V2)  components. 


3.2.  Moment  of  momentum  equation 


It  is  convenient  to  write  the  moment  of  momentum 
equation  [1] 

J'^r=Mb+Mh’  (8) 

in  the  M- coordinate  system  with  the  cylinder’s  angular 
velocity  components  (co l5co2,co3)  defined  by  (4).  Here, 
Mb  and  Mh  are  the  buoyancy  and  hydrodynamic  force 
torques.  The  buoyancy  force  induces  the  moment  in 
the  j M  direction  if  the  COM  does  not  coincide  with  the 
COV  (i.e.,  x  *  0), 

Mb  =  nxpwg  cos  (9) 


In  the  M-coordinate  system,  the  tensor  J  for  the  axially 
symmetric  cylinder  is  represented  by  a  diagonal  matrix 


J  = 


J  i  0  0 

0  J2  o 
0  0  J3 


(10) 


where  /2,  and  /3  are  the  moments  of  gyration  about 
the  center  of  mass.  The  gravity  force,  passing  the  center 
of  mass,  does  not  induce  the  moment. 


4.  Hydrodynamic  force  and  torque 

The  hydrodynamic  force  (Fh)  and  torque  (Mh)  are 
easily  calculated  in  the  F-coordinate  system  using  exist¬ 
ing  formulas.  After  calculation,  the  hydrodynamic  force 
(Fh)  should  be  transformed  from  the  F-coordinate  to  the 
F-coordinate  before  substituting  into  the  momentum 
Eq.  (7),  and  the  hydrodynamic  torque  (Mh)  should  be 
transformed  from  the  F-coordinate  to  the  M- coordinate 
before  substituting  into  the  moment  of  momentum 
equation  (8)  for  solutions  [3]. 

4.1.  Drag  force 

The  drag  force  consists  of  three  parts:  (a)  along  iF 
(along-cylinder)  drag-force  (Fdl),  (b)  along  (across- 
cylinder)  drag  force  (Fd2),  and  (3)  along  kp  drag 
force  (Fd3).  Let  (Cdi,Cd2)  be  the  drag  coefficients 
along-  and  across-cylinder  directions  (Reynolds  number 
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dependent).  The  drag  force  coefficients  are  calculated  on 
the  base  of  steady  flow,  it  is  different  from  the  fluid 
around  an  accelerated  solid  body.  The  added  mass  cor¬ 
rection  is  represented  by  the  ratios  (/i,  /2,  /3)  in  the  three 
directions  of  the  F-coordinate  system. 

The  drag  force  along-iF  is  calculated  by 

Fdi  =  \FCf\V\,  Cfi  =  Cdi  —  ^  lvil>  (H) 

where  Cd i  is  the  drag  coefficient  in  the  along-cylinder 
direction  and  less  dependent  on  the  axial  Reynolds  num¬ 
ber  (Re)  when  Re  >  104,  but  dependent  on  the  cylinder’s 
aspect  ratio,  <5  =  Lid.  An  empirical  formula  is  used  for 
calculating  Qi  [6], 

{1.0,  if  <5  >  8, 

0.75  +  <5/32.1934  +  0.09612/<52,  if  8  ^  <5  >  0.5, 
1.15,  if  <5  0.5. 

(12) 


Substitution  of  (4)  and  (5)  into  (11)  leads  to 


Fdi  —  —  Ctdi/jc 


u 

Uyj 

\ 

V 

- 

Vw 

w 

_Ww_ 

/ 

5  Lx  If  ip  5 


where  the  superscript  T’  denotes  the  transpose,  and  Ctd i 
is  the  drag  coefficient  along  the  cylinder  (in  the  ip  direc¬ 
tion)  for  the  relative  motion  between  COM  and  water 
(or  called  the  translational  drag). 

The  drag  force  along- is  calculated  by 


Fd?  = 


\\X \cd2(v2)2 

J-k-X  ^ 


-t -x  -  O+fh 

=  (Qd2^2  +/rd2)jF' 


df] 


Jf 


(13) 


where 


V'2(n)  =  V2- 

is  the  water-to-cylinder  velocity  at  the  surface  in  the  \F 
direction  and 


Ctd2 


Cd2T<7 


Pw 

6+ fi) 


(14a) 


Cd2  —  < 


1.9276  +  8 /Re, 
if  Re  <  12 
1.261  +  16/Re, 
if  12  <  Re  <  180 
0.855  +  89 /Re, 
if  180  <  Re  <  2000 
0.84  +  0.00003tfe, 
if  2000  <  Re  <  12,000 
1.2  —  4/S 

if  12, 000  <  Re  <  150, 000,  <5^10 
0.835  -  0.35/(5, 

if  12, 000  <  Re  <  150, 000, 2  <  <5  <  10 
0.7-0.08/5, 

if  12, 000  <  Re  <  150, 000,  S  <  2 
1.875  -  0.00000457?e, 
if  150,000  <7?e<  350,000 
1/(641550/^  +  1.5), 
if  Re  >  350, 000. 


Substitution  of  (4)  and  (5)  into  (13)  leads  to 


(15) 


Fd2  —  —Ctdlly 


u 

uw 

\ 

V 

- 

Vw 

w 

_Ww_ 

J 

+  /rddjh  '  lyy  ~  IfIf' 


(16) 

The  angular  velocity  (co'2)  around  jF  causes  non-uniform 
water-to-cylinder  velocity  in  the  kp  direction 

f3  =  o)'in.  (17) 

The  drag  force  along-kp  is  calculated  by 


Fd3  — 


^  ,  /  Pw 

2  dI  (iTTi) 


'Hi! 


rj2dri 


f . 

J+x 


rj2drj 


—  /rcD^F, 


(18) 


where 


is  the  translational  drag  coefficient  and 

fm  =  Cd2W(TTp)  &  + 


(14b) 


/rd3  =  -2cd2^^Z(3T2+4Z2)|a;'|(0',  (19) 

is  the  rotational  drag  force  in  the  kp  direction. 


is  the  rotational  drag  force  in  the  jF  direction.  Here, 
Cd2  is  the  drag  coefficient  for  the  across-cylinder  direc¬ 
tion.  An  empirical  formula  is  used  for  calculating  Cd2 
[7,6] 


4.2.  Lift  force 

The  water-to-cylinder  velocity  determines  the  lift 
force  [1] 
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F,  = 


F  [  v'lWidn 

1  J-h-i 


V  2 


(20) 


where  Q  is  the  lift  coefficient.  An  empirical  formula  is 
used  for  calculating  Q  [8], 

if  coi d/V2  ^  8, 


C 1  = 


=  t05Xd/V2, 

~  \  8  +  0 A2(co1d/V2  -  8),  if  cold/V2  >  8 


(21) 


and  Ct i  is  the  translational  lift  coefficient.  Substitution  of 
(4),  (5),  (14)  into  (20)  leads  to 


u 

uw 

\ 

V 

- 

Vw 

w 

_Ww_ 

J 

Fl  =  —Ct\Iy: 


where 
/rl  =  Ql 

is  the  rotational  lift  force. 


+  /rlkF,  Iyz  —  kFjp, 


(22) 


the  drag  force  causes  a  torque  on  the  cylinder  in  the 
kp  direction 


M, 


d3 


(h~x  1 
J+-*  2 


2C^(TTA) 


—  —  \Cm2(o'}  +  M3]  k/r, 

Pw 


( V2  —  a/3f/)2f/ffi7 
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l\2<- 


(1  +/r) 


1 
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M3^Cd2d(T^^X-  (29) 

The  lift  force  exerts  a  torque  on  the  cylinder  in  the  \F 
direction 


Mi,  = 


-/ 


1  0 

2Qd( TT« 


(F2  —  co'3o)ridri 


3f 


—  [Cmia)3  +  Mi]  jF, 


(30) 


4.3.  Hydrodynamic  torque 


For  an  axially  symmetric  cylinder,  the  hydrodynamic 
torque  in  the  iF  direction  is  not  caused  by  the  drag  and 
lift  forces,  but  by  the  viscous  fluid.  The  moment  of  the 
viscous  force  of  steady  flow  between  two  rotating  cylin¬ 
ders  with  the  common  axis  is  calculated  by  [5] 

2  2 

M  =  47l/|-y — (0J1  —  <X>o),  (23) 

ri  ro 

where  (r1?r0)  and  (coi,co0)  are  the  radii  and  angle  velo¬ 
cities  of  the  inner  and  outer  cylinders;  p  is  the  viscosity. 
Moment  of  the  viscous  force  on  one  rotating  cylinder  is 
the  limit  case  of  the  two  rotating  cylinders  as  r0  — >  oc, 
co0  =  0.  The  moment  of  the  viscous  force  around  iF  is 
calculated  by 

IVEcQ  CmiCOilp,  Cml  =  TCpLd  .  (^^) 

When  the  cylinder  rotates  around  )F  with  the  angular 
velocity  co'2 ,  the  drag  force  causes  a  torque  on  the  cylin¬ 
der  in  the  jF  direction 


Mh2  = 


,  I  ,  I  1 

-^kl  /  2 


2  Cad ( 1  ^ 


rj2\rj\dq 


3f 


—  [^2^2]  3f  5 


(3 m2  —  ^  C^d 


(h 


(l+f[)\l6L  +2LX 


(25) 

(26) 


where  /r  is  the  added  mass  factor  for  the  moment  of 
drag  and  lift  forces.  If  the  water-to-cylinder  velocity  or 
the  cylinder  mass  distribution  is  non-uniform  ( %  ^  0), 


c-sC,m(TT/»i 
"'^(TT T)lv^ 


(31) 


5.  Composition  of  model  cylinders 

5.1.  Description 

Three  model  cylinders  were  used  for  the  drop  exper¬ 
iment  at  the  Naval  Postgraduate  School  swimming  pool. 
All  have  the  same  diameter  of  0.04m,  but  different 
lengths  of  0.152,  0.121  and  0.091m.  The  bodies  are  con¬ 
structed  of  rigid  plastic  with  aluminum-capped  ends.  In¬ 
side  each  was  a  threaded  bolt,  running  lengthwise  across 
the  cylinder,  and  an  internal  mass  (Fig.  3).  The  internal 
cylindrical  mass  made  of  copper  is  used  to  vary  the  cyl¬ 
inder’s  center  of  mass  and  could  be  adjusted  fore  or  aft. 

The  model  cylinder  is  composed  of  six  uniform  cylin¬ 
drical  parts  (Fig.  4):  (1)  a  plastic  hollow  cylinder,  (2)  an 
aluminum-capped  left  end  solid  cylinder,  (3)  an  alumi¬ 
num-capped  right  end  solid  cylinder,  (4)  a  cylindrical 
threaded  rod,  (5)  a  cylindrical  threaded  bolt,  and  (6) 
an  adjustable  copper  cylindrical  mass  with  the  distance 
(2)  between  its  COMP  and  the  geometric  center  of  the 
model  cylinder  (B).  The  mass  in  each  cylindrical  part 
is  uniformly  distributed.  Therefore,  the  center  of  mass 
for  the  part  (COMP)  is  located  at  its  center  of  volume. 
The  geometric  characteristics  such  as  radius  (or  outer 
and  inner  radii),  length,  COMP  location  are  listed  in 
Tables  1  and  2.  The  density  of  each  part  is  listed  in  Table  3. 
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Fig.  3.  Internal  components  of  the  model  cylinder:  (1)  an  plastic 
hollow  cylinder;  (2)  an  aluminum-capped  left  end  solid  cylinder;  (3)  an 
aluminum-capped  right  end  solid  cylinder,  (4)  a  cylindrical  threaded 
rod;  (5)  a  cylindrical  threaded  bolt,  and  (6)  an  adjustable  copper 
cylindrical  mass. 
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Fig.  4.  Internal  structure  of  the  model  cylinder. 


Table  1 

Characteristics  of  cylindrical  parts 


Cylindrical  Mass 
parts 

Radius/outer 
and  inner  radii 

Length 

Location  of 

COMP  in  the 
M-coordinate 

1 

mi 

(Ru  Ri) 

L-2h 

1 

2 

m2 

Ri 

h 

LI 2  -  h/2  +  x 

3 

m3 

Ri 

h 

LI2  -  lxi 2  -  x 

4 

m4 

Ri 

L  —  2/] 

1 

5 

m5 

(■ R2 ,  R3) 

h 

Lt  2  -h-  m  - 1 

6 

m6 

(Ri,  R3) 

h 

^  +  z 

Table  2 


Geometric  parameters  of  the  cylindrical  parts  (unit:  10  2  m) 


Cylinder 

r2 

*3 

h 

h 

h 

L 

1 

2.0 

1.2 

0.3 

0.6 

1.9 

2.7 

15.2 

2 

2.0 

1.2 

0.3 

0.6 

1.9 

2.7 

12.1 

3 

2.0 

1.2 

0.3 

0.6 

1.9 

2.7 

9.1 

Table  3 


Density  (unit:  103  kgm  3  )  of  the  materials  for  the  model  cylinder 


Material 

Aluminum 

Copper 

Plastic 

Steel 

Density 

2.70 

8.93 

1.16 

7.80 

Table  4 

Moments  of  gyration  of  the  six  cylindrical  parts 


Cylindrical  part 

Ji 

Jl  =  J3 

1 

f{R\+R\) 

^R\  +  3R\  +  {L-2hf\ 

2 

rm  p2 

2  ^1 

®(3  R\  +  l\) 

3 

m  n2 

2  Kl 

1(3  R\  +  l\) 

4 

nt4  p  2 

^[3R23  +  (L-2h)2] 

5 

fRl 

f2(2R\  +  2R\  +  l\) 

6 

f(.Rl+Rl) 

^  (3ffl  +  3ff|  +  /|) 

5.2.  Moments  of  gyration 

The  moments  of  gyration  are  defined  by  von  Mises  [1] 


J  i  = 


/  {y1  +z2)d. 

Ji  =  J(x2  +y2)dm, 


M,  J  2 


/<-J +*=)<!'", 


(32) 


where  (x,y,z)  are  coordinates  in  the  M- coordinate  sys¬ 
tem.  Since  the  six  parts  of  the  cylinders  (or  hollow  cylin¬ 
ders)  have  uniform  mass  distribution,  the  moments  of 
gyration  for  the  part  (i)  are  easily  calculated  (Table  4). 
The  resultant  moments  of  gyration  (/1?  /2?  Jf)  for  the 
cylinder  are  computed  by 

j, = 


i=  1 

6 

J2=J3  =  +  miX2  +  m2 

i=  1 

( L-h 


L-h 


2  +  X J  +mAp  1-  l\- 


+  +  x)2- 

The  location  of  COM  (/)  is  determined  by 
[m5(L/2  -  l\  ~  h/2)  -  m62] 


1 


6 

E  mj 

i=  1 


(33) 

(34) 


which  indicates  how  the  adjustable  weight  determines 
the  location  of  COM  for  the  model  cylinder. 


6.  Similarity  for  experimental  design 

Most  Russian  sea  mines  usually  have  cylindrical 
geometry.  Small  and  large  aircraft-laid  mines  AMD-1- 
500  and  AMD- 1-1 000,  were  developed  in  1942.  Later, 
they  are  ranked  high  among  their  best  foreign  counter¬ 
parts.  At  the  end  of  WWII  (in  1945),  the  upgraded  ver¬ 
sions  of  these  mines  appeared.  They  were  designated 
AMD-2-500  and  AMD-2- 1000,  respectively.  These 
mines  are  designed  to  destroy  surface  warships  and 
other  vessels,  as  well  as  submarines.  Weight  and  size 
of  the  small  mine  are  similar  to  those  of  the  FAB-500 
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aerial  bomb  (500kg,  d  ~  0.45  m),  and  of  the  large  mine 
(1500kg,  d~  0.533m)  to  those  of  the  FAB-1500  aerial 
bomb.  Standard  aircraft  attachment  points  were  used 
to  suspend  the  mines.  Bottom  mines  (except  for  the 
AMD- 1-500)  can  be  scattered  from  surface  ships  and 
motor  torpedo  boats,  while  the  two  classes  of  the  large 
mines  can  be  deployed  from  the  above  platforms  and 
from  submarines.  The  aspect  ratio  (Lid)  is  always  much 
smaller  than  10  [9]. 

6.1.  Geometric  similarity 

Consider  the  shallow  water  scenario  of  naval  opera¬ 
tion  with  water  depths  of  12. 19-60. 96m  (i.e.,  40-200ft) 
[10].  Our  goal  was  to  choose  a  scale  that  was  somewhat 
representative  of  the  real  world  ratio  of  water  depth  to 
mine  length,  but  at  the  same  time  would  be  large  enough 
to  film  and  would  not  damage  the  pool’s  bottom.  The 
model  cylinders  were  based  on  the  realistic  assumption 
that  a  full  mine  with  length  of  3  m  is  laid  in  water  depths 
of  45m,  thus  producing  a  15:1  ratio.  The  depth  of  the 
pool  is  2.4m.  From  this  ratio,  the  length  ( L )  of  the 
model  mine  is  chosen  as  0.152  m.  The  addition  of  a 
0.121  and  0.091m  length  allowed  for  later  comparison 
of  the  sensitivity  of  water  phase  trajectory  to  the  ratio 
of  mine  length  over  diameter.  The  outer  radius  of  the 
model  mine  is  0.02  m.  The  aspect  (length/diameter)  ratios 
(Lid)  are  3.8,  3.025,  and  2.28. 

6.2.  Similarity  of  buoyancy  effect 

Usually,  the  density  ratio  (pi pw)  of  a  full  size  mine  is 
1.8.  The  model  cylinders  with  lengths  of  0.152,  0.121, 
0.091m  have  densities  of  (1.69, 1.67, 1.88)  x  103kgm-3, 
and  the  corresponding  density  ratios  of  (1.70,1.68, 
1.88).  In  each  of  the  three  model  cylinders,  the  location 
of  the  weight  (i.e.,  the  value  of  2)  is  adjustable.  Use  of 
(34)  location  of  the  COM  (^-value)  can  be  determined 
(Table  5).  The  positive  ^-value  indicates  that  COM  is 
below  COV,  and  the  negative  ^-value  indicates  that 
COM  is  above  COV.  The  maximum  ratio  of  yJL  is 
about  0.10,  0.08,  and  0.06  for  the  three  model  cylinders. 
Thus,  the  buoyancy  force  and  torque  are  similar 
between  the  model  cylinders  and  the  full  size  mines. 


6.3.  Hydrodynamic  similarity 

Hydrodynamic  force  and  torque  depend  on  the  drag 
coefficients  (Cd i,  Cd2)  and  lift  coefficient  Q.  The  drag 
coefficients  (12)  and  (15)  depends  on  the  aspect  ratio  S 
and  axial  Reynolds  number  (Re) 


where  v  =  0.18  x  10-5  m1 2  s-1,  is  molecular  viscosity  of 
the  water.  The  diameter  of  the  model  cylinder 
(d  =  0.04  m)  is  around  1/10  of  the  full  size  mine.  Usually, 
the  model  cylinders  (or  full  size  mines)  accelerate  after 
entering  the  water  surface.  The  axial  Reynolds  number 
increases  with  time.  For  the  shallow  water  scenario,  a 
full  size  mine  may  have  a  maximum  speed  just  few  times 
larger  than  that  of  the  model  cylinder.  For  F~2ms_1, 
the  axial  Reynolds  number  is  0.44  x  105  for  the  model 
cylinders  and  0.44  x  106  for  the  full  size  mine.  Since  it 
is  likely  Re>  104  for  the  model  cylinder  (or  the  full  size 
mine)  falling  through  the  water  column,  the  drag  coeffi¬ 
cient  along  the  cylinder  (Cdi)  depends  only  on  the  aspect 
ratio  5  [see  (12)],  and  therefore  is  comparable  for  the 
model  cylinder  and  for  the  full  size  mine  due  to  similar 
aspect  ratios. 

The  drag  coefficient  across  the  cylinder  (Cd2)  has  a 
more  complicated  relationship  than  Cdi  [see  (15)]. 
Although  the  maximum  axial  Reynolds  number  is  at 
least  10  time  larger  for  the  full  mines  than  for  the  model 
cylinders,  it  is  possible  the  axial  Reynolds  numbers  for 
the  full  size  mine  and  model  cylinder  fall  into  the  same 
interval  as  indicated  in  (15)  in  descending  through  the 
water  column.  The  phenomenon  observed  for  the  model 
cylinders  may  be  extrapolated  into  the  full  size  mines 
with  careful  consideration  of  the  range  of  the  axial 
Reynolds  number. 


7.  Cylinder  drop  experiment 

A  cylinder  drop  experiment  was  conducted  at  the 
NPS  swim  pool  in  June  2001.  The  purpose  of  the  experi¬ 
ment  is  to  collect  data  about  cylinder’s  motion  in  the 
water  column  for  various  combinations  of  the  model 


Table  5 

Physical  parameters  of  the  model  cylinders 

Cylinder  Mass  (kg)  L  (m)  Volume  (10-6m3)  pm  (103kgm-3)  J\  (10 4kgm2)  %  (10 2m)  J2  (J3)  (10 4kgm2) 

1  0.3225  0.152  191.01  1.69  0.3305  0.00  6.0879 

0.74  5.7830 

1.48  6.2338 

2  0.2542  0.121  152.05  1.67  0.2713  0.06  3.4246 

0.53  3.2065 

1.00  3.3126 

3  0.2153  0.0912  114.61  1.88  0.2350  0.00  1.6952 

0.29  1.5775 

0.58  1.5568 
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Fig.  5.  Experimental  equipments:  (A)  drop  angle  device,  (B)  cylinder 
injector,  (C)  infrared  light  sensor,  (D)  output  to  universal  counter,  and 
(E)  cylinders. 

mine  parameters.  The  purpose  is  to  detect  the  character¬ 
istics  of  the  cylinder’s  movement  in  the  water  column 
and  to  collect  the  data  of  the  location,  orientation. 

It  basically  consisted  of  dropping  each  of  three  model 
cylinders  into  the  water  where  each  drop  was  recorded 
underwater  from  two  viewpoints.  Fig.  5  depicts  the 
overall  setup.  The  controlled  parameters  for  each  drop 
were:  Lid  ratio,  /-value,  initial  speed  (F^),  and  drop 
angle.  The  ^-coordinate  system  is  chosen  with  the  origin 
at  the  corner  of  the  swimming  pool  with  the  two  sides  as 
x-  and  y-axis  and  the  vertical  z-axis.  The  initial  injection 
of  cylinders  was  in  the  (y,z)  plane  (Fig.  6). 

7.7.  Initial  speed 

Initial  speed  (F^)  was  calculated  by  using  the  volt¬ 
age  return  of  an  infrared  photo  detector  located  at  the 
base  of  the  cylinder  injector.  The  infrared  sensor  pro¬ 
duced  a  square  wave  pulse  when  no  light  was  detected 
due  to  blockage  caused  by  the  cylinder’s  passage.  The 


length  of  the  square  wave  pulse  was  converted  into  time 
by  using  a  universal  counter.  Dividing  the  cylinder’s 
length  by  the  universal  counter’s  time  yielded  F(m). 
The  cylinders  were  dropped  from  several  positions  with¬ 
in  the  injector  mechanism  in  order  to  produce  a  range  of 
The  method  used  to  determine  F^m)  required  that 
the  infrared  light  sensor  be  located  above  the  water’s 
surface.  This  distance  was  held  fixed  throughout  the 
experiment  at  10cm. 

7.2.  Drop  angle 

The  drop  angle  (initial  value  of  was  controlled 
using  the  drop  angle  device.  Five  screw  positions 
marked  the  15°,  30°,  45°,  60°,  and  75°.  The  drop  angles 
were  determined  from  the  lay  of  the  pool  walkway, 
which  was  assumed  to  be  parallel  to  the  water’s  surface. 
A  range  of  drop  angles  was  chosen  to  represent  the 
various  entry  angles  that  air  and  surface  laid  mines  exhi¬ 
bit  in  naval  operation.  This  range  produced  velocities 
whose  horizontal  and  vertical  components  varied  in 
magnitude.  This  allowed  for  comparison  of  cylinder  tra¬ 
jectory  sensitivity  with  the  varying  velocity  components. 

7.3.  Methodology 

For  each  drop  the  cylinder  was  set  to  a  /- value.  For 
positive  /- value,  the  cylinders  were  placed  into  the  injec¬ 
tor  so  that  the  COM  was  located  below  the  geometric 
center.  For  negative  /-value,  the  COM  was  located 
above  the  geometric  center  to  release.  A  series  of  drops 
were  then  conducted  in  order  of  decreasing  mine  length 
for  each  angle.  Table  6  indicates  number  of  drops  con¬ 
ducted  for  different  drop  angles  and  /-value  for 
Ltd  =3.85.  Number  of  drops  for  other  Ltd  ratios 
(3.025,2.28)  is  comparable  to  that  for  Ltd  ratio  of 
3.85.  All  together  there  were  230  drops.  Each  video  cam¬ 
era  had  a  film  time  of  approximately  one  hour.  At  the 
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Fig.  6.  Top  view  of  the  cylinder  drop  experiment.  The  two  video  cameras  are  in  the  middle  of  the  water  column  (water  depth  2.4m). 
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Table  6 

Number  of  drops  conducted  for  different  drop  angles  and  ^-values  for 
Lid  =  3.85 


<4to) 

15° 

30° 

45° 

o 

O 

75° 

12 

13 

15 

15 

15 

12 

1\ 

9 

15 

15 

15 

9 

lo 

12 

14 

15 

18 

6 

l-l 

0 

6 

6 

6 

0 

1-2 

2 

6 

6 

0 

0 

end  of  the  day,  the  tapes  were  replayed  in  order  to  deter¬ 
mine  clarity  and  optimum  camera  position. 

8.  Data  retrieval  and  analysis 

8.1.  Data  retrieval 

Upon  completion  of  the  drop  phase,  the  video  from 
each  camera  was  converted  to  digital  format.  The  digital 
video  for  each  view  was  then  analyzed  frame  by  frame 
(30  Hz)  in  order  to  determine  the  mine’s  position  in  the 
(pc,z)  and  (y,z)  planes.  The  mine’s  top  and  bottom 
positions  were  input  into  a  MATLAB  generated  grid, 
similar  to  the  ones  within  the  pool.  The  first  point  to  im¬ 
pact  the  water  was  always  plotted  first.  This  facilitated 
tracking  of  the  initial  entry  point  throughout  the  water 
column.  The  cameras  were  not  time  synchronized;  thus, 
the  first  recorded  position  corresponded  to  when  the  full 
length  of  the  mine  was  in  view. 

8.2.  Source  of  errors 

There  were  several  sources  of  error  that  hindered  the 
determination  of  the  cylinder’s  exact  position  within  the 
water  column.  Locations  above  or  below  the  camera’s 
focal  point  were  subjected  to  parallax  distortion.  Placing 
the  cameras  as  far  back  as  possible,  while  still  being  able 
to  resolve  the  individual  grid  squares,  minimized  this 
error.  Second,  the  background  grids  were  located  behind 
the  cylinder’s  trajectory  plane.  This  resulted  in  the  cylin¬ 
der  appearing  larger  than  normal.  This  error  was  mini¬ 
mized  by  neglecting  that  data  if  the  plotted  points 
exceed  the  particular  cylinder’s  length.  Third,  an  object 
injected  into  the  water  will  generate  an  air  cavity.  This 
air  cavity  can  greatly  affect  the  initial  motion,  particu¬ 
larly  at  very  high  speeds  (hydro  ballistics).  The  air  cavity 
effect  was  deemed  to  be  minimal  due  to  the  low  inject 
velocities  used. 

8.3.  Data  analysis 

The  data  provided  by  each  camera  was  first  used  to 
produce  raw  two  2-D  plots  of  the  cylinder’s  trajectory. 
Next,  2-D  data  from  both  cameras  was  then  fused  to 


produce  a  3-D  history.  This  3-D  history  was  then  made 
non-dimensional  in  order  to  generalize  the  results.  The 
non-dimensional  data  was  used  to  generate  impact  scat¬ 
ter  plots  and  was  also  used  in  multiple  linear  regression 
calculations. 


9.  Experimental  results 

9.1.  Trajectory  patterns 

After  analyzing  the  3-D  data  set,  seven  trajectory  pat¬ 
terns  were  found.  The  plots  on  the  (y,z)  plane  were  cho¬ 
sen  for  trajectory  analysis,  as  this  plane  was  parallel  to 
the  direction  of  the  mine  drop.  Observed  trajectories 
(Fig.  7)  were  found  to  be  most  sensitive  to  ^-value,  drop 
angle  and  Ltd  ratio  (Table  7).  Dependence  of  the  trajec¬ 
tory  patterns  on  the  cylinders’  physical  parameters  and 
release  conditions  are  illustrated  in  Table  8  for  %  >  0 
and  Table  9  for  x  <  0. 

For  positive  ^-values  (nose-down),  as  the  distance  be¬ 
tween  COM  and  COV  increases  (i.e.,  x  increases),  the 
cylinder’s  trajectory  tends  to  follow  the  straight  pattern 
(Table  8).  As  the  distance  between  COM  and  COV  de¬ 
creases  (i.e.,  x  decreases),  the  cylinder’s  trajectory  tends 
towards  being  more  parallel  with  the  pool’s  bottom. 
At  steep  drop  angles,  the  cylinder  experiences  little  lat¬ 
eral  movement  and  tends  towards  a  straight  pattern. 
Additionally,  as  Lid  ratio  decreased  more  complex  tra¬ 
jectory  patterns  developed.  This  included  significant 
oscillation  about  the  vertical  axis  and  increased  lateral 
movement. 

For  negative  ^-values  (nose-up),  the  cylinder  will  flip 
only  once  and  almost  all  the  flips  occur  right  after  enter¬ 
ing  the  water  surface.  This  is  true  for  all  30  nose-up  drop 
( X  <  0)  cases  (Table  9).  As  the  distance  between  COM 
and  COV  increases  (i.e.,  \x\  increases),  the  cylinder’s  tra¬ 
jectory  tends  to  follow  the  flip-straight  pattern  (Table 
9).  As  the  distance  between  COM  and  COV  decreases 
(i.e.,  x  decreases),  the  cylinder’s  trajectory  tends  towards 
the  flip-straight-spiral  or  flip-spiral  pattern.  Similar  to 
the  nose-down  case,  as  Lid  ratio  decreased  more  com¬ 
plex  trajectory  patterns  developed. 

9.2.  Impact  attitude 

The  angle,  i fa  +  n/2,  is  the  impact  attitude  at  the  bot¬ 
tom  of  the  water  column.  The  mine  burial  is  largely 
determined  from  the  impact  attitude  of  the  cylinder. 
Cylinders  whose  impact  attitudes  are  perpendicular 
(x/j 2  =  90°)  to  the  sediment  interface  will  experience  the 
largest  degree  of  impact  burial  [1 1].  It  is  therefore  impor¬ 
tant  to  analyze  the  relationship  between  impact  attitude 
and  the  controlled  parameters,  drop  angle,  V^m\  Lid , 
and  x-  The  experiment  shows  that  both  Ltd  and  V^in) 
had  little  influence  on  impact  attitude.  The  impact 
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Flip 
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Flip-Flat 

L-  9cm;  Vi-  2  294m/s 


Y  (m) 


Y  (m) 


Fig.  7.  Cylinders’  track  patterns  observed  during  the  experiment. 


attitude  largely  depends  on  the  value.  The  histogram 
of  the  impact  angle  (Fig.  8)  shows  five  peaks  centered 
near  0°,  40°,  90°,  140°,  and  180°  corresponding  to  the 
COM  positions  (X-2,X-uXo,XuX2)  [he.,  (COM- 
2,  COM- 1 ,  COMO,  COM  1 ,  COM2)] . 

Although  drop  angle  is  not  the  most  influential 
parameter,  variations  did  induce  changes  in  impact  ori¬ 
entation.  As  drop  angle  increases,  the  likelihood  of  any 
lateral  movement  decreases.  This  allowed  for  impact  an¬ 
gles  that  are  more  vertically  orientated.  This  is  primarily 


due  to  the  fact  that  the  vertical  components  of  velocity 
are  greater  than  those  at  shallow  angles.  Thus,  the  time 
to  bottom  and  time  for  trajectory  alteration  is  less. 

10.  Numerical  simulation 

A  numerical  model  on  the  base  of  momentum  bal¬ 
ance  (7)  and  moment  of  momentum  balance  (8)  has  been 
developed  to  predict  the  cylinder’s  translation  velocity 
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Table  7 

Characteristics  of  trajectory  patterns 
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Trajectory  pattern 

Description 

Straight 

Cylinder  exhibited  little  angular  change 
about  z-axis.  The  attitude  remained  nearly 
parallel  with  z-axis  (±15°) 

Slant 

Cylinder  exhibited  little  angular  change 
about  z-axis.  The  attitude  was  45°  off 
z-axis  (±15°) 

Spiral 

Cylinder  experienced  rotation  about  z-axis 
throughout  the  water  column 

Flip 

Initial  water  entry  point  rotated  at  least  180° 

Flat 

Cylinder’s  angle  with  vertical  near  90°  for  most 
of  the  trajectory 

Seesaw 

Similar  to  the  flat  pattern  except  that 
cylinder’s  angle  with  vertical  would  oscillate 
between  greater  (less)  than  90°  and  less 
(greater)  than  90° — like  a  seesaw 

Combination 

Complex  trajectory  where  cylinder  exhibited 
several  of  the  above  patterns 

and  orientation  using  the  triple  coordinate  transform  [3]. 
The  numerical  model  predicts  the  motion  of  cylinder  in¬ 
side  the  water  column  reasonably  well.  Two  examples 
are  listed  for  illustration. 

Positive  x  (nose-down) .  Cylinder  #1  (X  =  0.152m, 
p  =  1.69  x  103kgm-3)  with  x  ~  0.0074  m  is  injected  to 
the  water  with  the  drop  angle  45°.  The  physical  para¬ 
meters  of  this  cylinder  are  given  by 

m  =  0.3225kg,  J\  =  0.3305  x  10-4kgm2, 

J2=J3  =  5.783  x  10“4kgm2.  (36a) 


Impact  Angle  Frequency  of  Occurrence 


Undersea  cameras  measure  the  initial  conditions 

x0  =  0,  y0  =  0,  zq  =  0,  u0  =  0, 

vq  =  —  1.55ms-1,  wo  =  —  2.52ms-1, 

*Aio  =  0,  4^20  =  60  ,  i// 30  =  —95  , 

ojio  =  0,  oj2()  =  0.49  s-1,  co30  =  0.29  s_1.  (36b) 

Substitution  of  the  model  parameters  (36a)  and  the  ini¬ 
tial  conditions  (36b)  into  the  numerical  model  (see  [3]) 
leads  to  the  prediction  of  the  cylinder’s  translation  and 
orientation  that  are  compared  with  the  data  collected 
during  the  experiment  at  time  steps  (Fig.  9).  Both  model 
simulated  and  observed  tracks  show  a  slant-straight 
pattern. 


Table  8 

Dependence  of  trajectory  patterns  on  input  conditions  for  nose-down  dropping  (/  >  0) 


Cylinder  length  (m) 

0.152 

0.121 

0.0912 

1  (m) 

0.0148 

0.01 

0.0058 

Drop  angle  15° 

Straight  (1)  slant-straight*  (3) 

Straight  (1),  spiral  (1)  slant-straight*  (2) 

Spiral*  (2)  straight-slant  (1) 
slant-straight  (1) 

Drop  angle  30° 

Straight  (1)  slant-straight*  (4) 

Slant  (1),  spiral  (1)  straight  (1) 
slant-straight*  (2) 

Spiral*  (5) 

Drop  angle  45° 

Slant*  (2),  straight  (1)  slant-straight 
(1)  straight-spiral  (1) 

Straight  (1)  spiral  (2)  straight-spiral 
(1)  slant-straight  (1) 

Spiral*  (4)  slant-spiral  (1) 

Drop  angle  60° 

Straight**  (5) 

Straight*  (3)  Straight-spiral 
(1)  straight-slant  (1) 

Spiral  (4)  straight-spiral  (1) 

Drop  angle  75° 

Straight**  (5) 

Straight  (2)  straight-spiral  (3) 

Spiral  (2),  slant  (1)  straight-spiral  (2) 

Here,  the  symbols  and  ***’  mean  the  dominant  and  only  patterns  in  the  given  bin. 


Table  9 

Dependence  of  trajectory  patterns  on  input  conditions  for  nose-down  dropping  (%  <  0) 


Cylinder  length  (m) 

0.1520 

0.1210 

0.0912 

i  (m) 

0.0148 

0.01 

0.0058 

Drop  angle  30° 

Flip-straight*  (3)  flip-slant  (1) 

Flip-straight  (2)  flip-slant  (2) 

Flip-straight-spiral*  (3)  straight-flip-seesaw  (1) 

Drop  angle  45° 

Flip-straight*  (2)  Flip-slant 

Flip-straight  (1)  flip-slant* 

Flip-straight  (1)  flip-spiral* 

(1)  flip-spiral-slant  (1) 

(2)  flip-straight-spiral  (1) 

(2)  flip-straight-spiral  (1) 

Drop  angle  60° 

Flip-straight  (1)  straight-flip  (1) 

Flip-straight  (1)  slant-flip-slant  (1) 

Flip-spiral-seesaw  (1)  flip-spiral  (1) 

Here,  the  symbol  **’  means  the  dominant  pattern  in  the  given  bin. 
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Experiment  Model 


Fig.  9.  Movement  of  cylinder  #1  (L  =  0.152m,  p  =  1.69  x  103kgm-3) 
with  x  =  0.0074m  and  drop  angle  45°  obtained  from:  (a)  experiment 
and  (b)  recursive  model. 


Negative  /  (nose-up).  Cylinder  #2  ( L  -  0.121m,  p  - 

I. 67  x  103kgm-3)  with  y  =  —  0.01m  is  injected  to  the 
water  with  the  drop  angle  30°.  The  physical  parameters 
of  this  cylinder  are  given  by 

m  =  0.2542kg,  Jx  =  0.2713  x  10-4kgm2, 

J2=J3  =  3.313  x  10-4kgm2.  (37a) 

Undersea  cameras  measure  the  initial  conditions 
x0  =  0,  y0  =  0,  zq  =  0,  u0  =  0, 
v()  =  —0.75m  s'1 ,  wo  =  —0.67 m s1 , 

*Aio  =  0?  ^20  =  24  ^  ^30  =  ~96  j 

co i0  =  0,  co2o  =  —  5.08  s-1,  m30  =  0.15  s-1.  (37b) 

The  predicted  cylinder’s  translation  and  orientation  are 
compared  with  the  data  collected  during  the  experiment 
at  time  steps  (Fig.  10).  The  simulated  and  observed 
tracks  show  a  flip-spiral  pattern.  The  flip  occurs  at 
0.11s  (0.13  s)  after  the  cylinder  enters  the  water  in  the 
experiment  (model).  After  the  flip,  the  cylinder  spirals 
down  to  the  bottom. 

II.  Statistical  analysis 

Non-dimensional  velocity  data  are  got  using 

V  =  VgZv  *  .  (38) 

Multivariate  linear  regression  is  conducted  on  the  exper¬ 
imental  data  to  establish  relationships  between  the  phys¬ 
ical  parameters  of  the  cylinder  {Lid,  y),  initial  condition 
(F*H,  i/4m))  and  the  output  data  (temporally  varying) 


Experiment  Model 


Fig.  10.  Movement  of  cylinder  #2  (L  =  0.121  m,  p  =  1.67  x  103kgm-3) 
with  /  =  -  0.01m  and  drop  angle  30°  obtained  from:  (a)  experiment 
and  (b)  recursive  model. 


such  velocity  (w*,p*,w*)  and  attitude  i/A2)  at  time  t. 
Let  Y{t )  represent  the  output  variables.  The  regression 
equation  is  given  by 

Y{t)  =  p0(t)  +  faM™  +  P2(t)L/d  +  &  +  P4(t)X, 

(39) 

which  is  statistically  significant  on  0.01  confidence  level 
after  the  correlation  coefficient  test.  Fig.  11  shows  the 
temporally  varying  regression  coefficients  for  the  varia¬ 
ble  u*,  v*,  w*,  and  i//2-  For  the  cylinder’s  attitude  (i J/2) 
the  coefficient  f)4{t)  is  much  larger  than  the  other  coeffi¬ 
cients,  which  indicates  that  the  cylinder’s  orientation 
(versus  the  vertical  direction)  is  mainly  determined  by 
the  location  of  COM  (i.e.,  the  value  of  y). 

To  show  the  dependence  of  impact  of  the  cylinder  at 
the  bottom  (i.e.,  horizontal  location  relative  to  the  cylin¬ 
der’s  surface  entry  point,  orientation,  velocity),  the  mul¬ 
tivariate  regression  between  the  input  non-dimensional 
parameters:  il/2n\  Lid ,  F*(m),  and  %,  and  the  final  state 
(i.e.,  impact  on  the  bottom)  variables  such  as  the  hori¬ 
zontal  position  of  COM  (xm,ym),  the  velocity  of  COM 
(w*,p*,w*)  and  the  attitude  i j/2.  Let  Z  represent  the  out¬ 
put  variables.  The  regression  equation  is  given  by 

Z  =  oc0  +  aii/4m)  +  oc2L/d  +  a3F*(m)  +  a4^,  (40) 

with  the  regression  coefficients  (a0,ai,a2,  a3,a4)  for  the 
output  parameters  given  by  Table  10.  The  impact  alti¬ 
tude  at  the  bottom  is  given  by 

t/4bot)  =  103  -  13.4t/4in)  -  0.501  L/d  +  1.045r(in)  +  472/, 
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Fig.  11.  Temporally  varying  regression  coefficients. 


Table  10 


Regression  coefficients  of  predicting  bottom  impact 


^2 

u* 

v* 

w* 

a0 

(4in)) 

84.22 

-0.4675 

-0.1074 

-1.6268 

a  i 

0.354 

0.002 

-0.0002 

0.0006 

a  2 

(Lid) 

-3.329 

0.216 

0.0581 

0.2251 

a3 

(V*(in)) 

-1.457 

-0.0063 

-0.0331 

0.0615 

a4 

(x) 

625.65 

0.591 

-0.0861 

0.1657 

where  \\j2  is  in  degree.  The  regression  coefficients  are  1- 
3  orders  of  magnitude  smaller  for  (Lid)  and  F*(m)  than 
for  \j/^  and  x>  which  indicates  that  Lid  and  initial 
water  entry  speed  have  little  influence  on  cylinder’s  im¬ 
pact  attitude  on  the  bottom.  For  %  =  0,  the  cylinder  is 
almost  parallel  to  the  bottom.  For  x~2  and  X2  cases 
(large  \x\  value),  the  cylinder  is  almost  vertical  to  the 
bottom. 


12.  Conclusions 

(1)  Movement  of  a  falling  cylinder  in  water  column  is 
a  highly  nonlinear  process.  Six  trajectory  patterns 
(straight,  spiral,  flip,  flat,  seesaw,  combination) 
are  detected  from  the  experiment.  The  transi¬ 
tion  between  patterns  depends  on  the  initial  condi¬ 
tions  (drop  angle  xj/™)  and  initial  speed  V^in) 
and  the  physical  parameters  of  the  cylinder  (such 
as  Lid  ratio,  ^-value).  All  the  nose-up  drops 
(X  <  0,  30  drops)  shows  that  the  cylinder  flips  only 
once  once  for  negative  ^-values.  The  flat  pattern 
occurs  usually  for  x  =  0  (i.e.,  COM  coincides  with 
COV). 

(2)  The  experiment  shows  that  both  Lid  and  V(m)  had  lit¬ 
tle  influence  on  cylinder’s  impact  attitude  on  the  bot¬ 
tom.  The  x-value  is  most  important  to  determine  the 
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impact  attitude.  For  %  =  0,  the  cylinder  is  almost  par¬ 
allel  to  the  bottom.  For  x~i  and  Xi  cases,  the  cylinder 
is  almost  vertical  to  the  bottom. 

(3)  The  experiment  provides  a  data  set  of  temporally  var¬ 
ying  cylinder’s  orientation,  COM  position  and  velo¬ 
city  for  various  input  conditions  such  as  cylinder’s 
physical  parameters  and  initial  conditions.  The  data 
are  useful  for  model  development  and  validation. 

(4)  The  cylinder’s  movement  during  the  experiment  is 
well  simulated  using  a  numerical  model  based  on 
the  momentum  and  moment  of  momentum  balances 
with  triple  coordinate  transform. 
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